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New vistas of the meson structure in QCD from low to high energies 
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This talk presents issues pertaining to the quark structure of the pion within QCD, both from the theoretical 
and from the experimental point of view. We review and discuss the pion-photon transition form factor and the 
pion's electromagnetic form factor vs. corresponding experimental data from the CLEG Collaboration and the 
JLab. We also examine the extent to which recent high-precision lattice computations of the second moment 
of the pion's distribution amplitude conform with theoretical models. Finally, we include predictions for the 
azimuthal asymmetry of the distribution in the polarized /i-pair-induced DY production employing various 
pion distribution amplitudes. 



1. Introduction 

Understanding the quark structure of the hght- 
est meson, the pion, — the prototype for a meson 
bound state — is arguably one of the most basic, 
albeit challenging, questions QCD is still facing 
even after decades of intense investigations. Inte- 
grating over transverse momenta up to some reso- 
lution scale ^Q, one gets from the pion's light-cone 
wave function the pion's leading twist-2 distribu- 
tion amplitude (DA), iPt^(x,^q), defined in terms 
of a nonlocal axial current: 

(0| d(z)7^75C(z, 0)u(0)K(P))U2^o = 

tUP^^f dxe"(-^Vr"'(a;,Mo), (1) 
Jo 

where x {x = 1 ~ x) is the longitudinal momen- 
tum fraction carried by the valence quark (an- 
tiquark) in the pion and the path-ordered ex- 
ponential C(z,0) = Ve^pl-ig dyH'^A'^^y)] 
ensures gauge invariance. Appealing to its 
renormalization-group properties [112] , we can ex- 
pand ip'^'"~^{x, Hq) in terms of its one-loop eigen- 
functions, alias the Gegenbauer polynomials, to 



obtain 



(2) 



*TaIk presented at International Workshop on e+e col- 
lisions from Phi to Psi, Frascati, Italy, 7-10 April 2008. 



Here ip'^^^x) = 6xx is the asymptotic pion DA 
and by virtue of the leptonic decay tt — s- ^i^v^ 
one has the normalization dx ip^""~'^{x, (Iq) = 

1. Relying only upon the first two Gegenbauer 
coefiicients, Eq. ([2]) can yield distinct profiles, as 
shown in Fig. [TJ We will see in the next section 
how two-photon processes can be used to resolve 
pion's dilemma and reveal its parton substructure 
in agreement with the experimental data and the 
latest lattice calculations. 

2. Pion DA, CLEO data, and lattice esti- 
mates 

The entire nonperturbative content of the pion 
DA is encoded in the expansion coefficients 
o-n (/^^ ) which in turn can be derived from the mo- 
ments 

(e^)^EE f dx{2x~l)''^,{x), (3) 
Jo 

where ^ = 2.t — 1, that decrease with increas- 
ing polynomial order to 0: (^^)^ ^ [3/(A^ + 
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Fi gure 1 . Pion's dilemma: to B (Bactrian camel) — 
dashed line (CZ) [3] and BMS "bunch" g]— or to D 
(Dromedary camel) — dotted line (asymptotic DA)? All 
curves are shown at the normalization scale /^^ = 1 GeV^ . 

1)(A^ + 3)]. The evolution behavior of a„(/.t^) is 
controlled by the ERBL equation [112] (for a ped- 
agogical exposition, see [H])- In our approach U 
we have determined a„ at a normalization scale 
= 1.35 GeV^ with the help of QCD sum rules 
with nonlocal condensates |6l7j . There are alter- 
native approaches, based, for instance, on local 
QCD sum rules [5], instantons |8l9ll0j . etc. 

The predictive power of these theoretical ap- 
proaches was challenged by the high-precision 
data of the CLEO Collaboration on the pion- 
photon transition in which one photon has 
a large virtuality Q^, while the other is nearly 
on shell. Kroll and Raulfs [T^ were the first to 
show that the popular CZ pion DA was over- 
shooting these data considerably. Other analy- 
ses, having recourse to light-cone sum rules with 
a spectral density obtained in the standard fac- 
torization scheme of perturbative QCD in leading 
order (LO) [13] and next-to-leading-order (NLO) 
|14I15I16I17] . followed, which established the fol- 
lowing facts (consult Fig. [2l drawn at the main 
scale, A*cLEO = (2-4 GeV) , probed in the CLEO 
experiment): (i) The CZ pion DA is outside the 
4ct error eUipse of the CLEO data, (ii) The 
asymptotic pion DA is outside the 3(t ellipse, 
while (iii) other proposed models close to that 
are at least 2a's off. (iv) The "bunch" of pion 
DAs derived from nonlocal QCD sum rules mostly 
overlaps with the la error ellipse, with the BMS 
model DA being entirely inside. 

The most recent lattice calculations [18] (larger 
band bounded by dashed lines) and [19] (narrower 
band within solid lines) in Fig. [5] support and en- 
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Figure 2. Comparison at Mcleo various pion DA 
models with the CLEO data in terms of the la ellipse 
(solid line) and recent lattice simulations, denoted by ver- 
tical dashed lines 1181 and solid ones 1191 . The symbols 
mark the models discussed in Tables [T] and |2] The slanted 
shaded rectangle represents the BMS DA "bunch" [4]. The 
dashed la ellipse corresponds to the inclusion of the twist- 
4 contribution to the pion DA via renormalons 1171 . The 
red bullet denotes the prediction derived from the holo- 
graphic model based on AdS/CFT duality (see last sec- 
tion). 
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Figure 3. Range of values of (5^)7r (Mcleq) (slanted bro- 
ken lines) intersecting with the la error ellipse (inner solid 
line) of the CLEO data and the latest lattice results on a2 
of the UKQCD-RBC Collaboration [19] (vertical shaded 
band) . 

hance these findings as regards the range of val- 
ues of the first Gegenbauer coefhcient a2. First, 
as one sees from this figure, they rather disfa- 
vor a relatively large twist-4 contribution to the 
pion DA, estimated with the help of renormalons 
[TT] (dashed la error eUipse). Second, the latest 
calculation [19| . with even smaller uncertainties 
than [18] , indicates a trend further away from the 
asymptotic DA, but still in compliance with the 
BMS results. 

Tables [1] and O show the values of the second 
and the fourth moment of various pion DAs at 
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the lattice reference scale /i^at = 4 GeV^. 
Table 1 



Predictions for {S,'^)-^ from various models (meth- 
ods) after NLO evolution from their proprietary 
normalization point to the scale //^at = ^ GeV^ . 



Source 


Min 


Mid 


Max 


UKQCD/RBC U2i 


0.252 


0.278 


0.304 


QCDSF/UKQCD [ig 


0.230 


0.269 


0.308 


m 


0.232 


0.285 


0.338 


[TSITB] (iCT CLEO) + 


0.240 


0.280 


0.320 




0.233 


0.248 


0.264 


m* 


0.233 


0.269 


0.305 


uni^ 




0.223 




m ^ 




0.212 




g] ^ 




0.211 








0.259 




m ^ 




0.229 




[23] T 




0.312 




m ■ 




0.344 




Asy ♦ 




0.199 




AdS/QCD • 




0.237 





In Fig. [31 we determine the range of values of 
the fourth moment (^^}7r (02,04!^^) (denoted by 
the slanted broken lines) that simultaneously ful- 
fill the CLEO data and the lattice constraints on 
{^'^)tt from [19]. For each error ellipse, there is 
some overlap with the vertical a2-band estimated 
in |19| . For the maximum and minimum of (^^)7r 
with respect to the Icr error ellipse, we find 

0.095 < (e*). (mL) < 0.134, (4) 

corresponding to (a^'" = 0.149, af"^ = -0.241) 
and (af^^^ = 0.307, af^^ = -0.205) . [Recall that 
^^4^asy _ 3/35.] It happens that these values al- 
most coincide with the intersection points (open 
circles) of the la error ellipse and the boundaries 
of the vertical band of the a2 lattice constraints. 

Figure [4| displays various theoretical predic- 
tions for (5^F-y.^_+^((3^) vs. in compari- 
son with the CELLO (diamonds, [53]) and the 
CLEO (triangles, [11]) experimental data, evalu- 
ated with the twist-4 parameter value (5xw-4 ^ 
0.19 GeV^ |15I16I17| . The other curves shown 



Table 2 



Predictions for (^'*)7r from various models (meth- 
ods) after NLO evolution from their proprietary 
normalization point to the scale ^^at = ^ GeV^ . 



Source 


Min 


Mid 


Max 


[15|16j (Icr CLEO) 


0.066 


0.115 


0.162 


01 


0.105 


0.108 


0.113 


m 


0.082 


0.116 


0.151 


[To' 




0.10 




m 




0.091 









0.094 




[2r 




0.121 




m 




0.104 




m 




0.178 




m 




0.181 




Asy 




0.085 




AdS/QCD 




0.114 





correspond to selected pion DAs: the asymptotic 
DA ipas (lower dashed line), ipcz (upper dashed 
line) [J, and two instanton-based models, viz., 
[8] (dotted line) and [10] (dash-dotted line). An 
important observation from this figure is that 
the shaded strip, which corresponds to the BMS 
"bunch" [1], becomes narrower at lower scales 
around 1 GeV^. The reason is that at such low 
scales, the form factor is dominated by its twist-4 
contribution, while the leading twist-2 part dies 
out. Moreover, Fig. |4| makes it clear that the 
low-Q^ CELLO data [23] exclude (p^''^ and clones, 
while the high-Q^ CLEO data [11] rule out ip^^. 
Figure [5] gives an illustration of the partial con- 
tributions to Q^F^ '^'^ , originating from different 
sources at /icLEO- ^ comprehensive account of 
these effects can be found in [T7] . 

3. Electromagnetic pion form factor 

The analysis of the pion's electromagnetic form 
factor involves on the nonperturbative side the 
BMS "bunch" of pion DAs in comparison with 
1^9^'''' and v^J^. On the perturbative side, a theo- 
retical scheme is used, which consists of expand- 
ing the form factor in terms of analytic images 
of the strong running coupling and its powers up 
to the NLO. The basis of this approach develops 
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Figure 4. Predictions for Q2^^,^_^^(q2) q2 
comparison with available experimental data (see text). 




Figure 5. Contributions to Q^FT*^"^ at the typical scale 
fi^ = 5.76 GeV^ of the CLEO data [U]. The next-to- 
next-to-loading order (NNLO) estimate is based on |25) . 
whereas the uncertainty owing to the (Ros)onance model 
used, was discussed in 1131 . 

from |26| and can account for more than one hard 
scale by incorporating into the "analytization" 
procedure all terms that contribute to the spec- 
tral density of the amplitude [57] . The bedrock of 
the approach was developed in |27I28I29I30I31I32| 
and its application to F^^'~'{Q'^) was considered 
in [33|34| (see also [35]). The results are displayed 
in Fig. [6] in comparison with experimental data. 
Note that the quantity shown comprises a soft 
non-factorizing contribution which dominates at 
currently accessible momentum transfers [33] . 

The main characteristics of this approach are 
(i) a strongly reduced sensitivity on the renormal- 
ization and the factorization scale, (ii) an undi- 
minished quality of precision in adopting different 
choices of renormalization schemes and scale set- 
tings, virtually eliminating the dependence varia- 
tions from scheme to scheme and scale to scale 
[34] . (iii) Another important finding is that, 
within such an analytic approach, the form-factor 
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Fi gure 6. Predictions for the scaled pion form factor cal- 
culated with the BMS "bunch" (shaded strip) |4] in NLO 
QCD analytic perturbation theory I33| . The dashed lines 
inside the strip restrict the area of predictions accessible 
to the asymptotic pion DA. The experimental data are 
taken from 1361 (diamonds) and 1371 (triangles). 



predictions (shaded strip in Fig. [5]) turn out to 
be very close to that computed with (fi'^^^ , albeit 
the underlying pion DA profiles are very different. 
This proved that what really matters is the be- 
havior of the pion DA at the endpoints a; — > 0, 1. 



4. Conclusions 

We have shown that the CLEO data on F''''*'' 
pose a veto to a variety of proposed models for 
the pion DA (see Fig. ^ and favor an endpoint- 
suppressed "B" -shaped pion DA — a Bactrian 
"camelino", the endpoint suppression being pro- 
vided by the vacuum nonlocality = 0.4 GeV^. 
The crucial 'missing link' is the determination of 
on the lattice. A value within the range 
[0.095 -0.134] (2 GeV), extracted from the CLEO 
data in conjunction with the lattice constraints 
on (^^),r from [THj (cf. Fig. [3]), would vahdate the 
claim that the pion DA is BMS-like. 

The Drell-Yan process 7t~N fi^fi^X for 
lepton-pair production with a large invariant 
mass provides an additional useful tool to 
probe and test different pion DAs in terms of 
azimuthal asymmetries, as Fig. [7] illustrates for 
the case of the kinematic variable /i. Overall, 
a rather good agreement of the BMS "bunch", 
derived from nonlocal QCD sum rules [?], with 
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Fi gure 7. Results for the angular distribution parameter 
^ as a function of Xu = Xr: for p = Qt/Q = 0.5. The 
shaded strip contains the results for the BMS "bunch" of 
DAs [4]. The solid line corresponds to the BMS model, 
the dotted solid denotes the result for , and the (red) 
dashed line is the prediction for the endpoint-dominated 
CZ DA [3]. One-loop evolution of the pion DAs to each 
measured value is included (data taken from [38]). 

available data was found [39]; though the exist- 
ing data cannot single out a particular model. In 
this respect, the planned COMPASS experiment 
may be of significant relevance. 

5. Section added: Holographic dual of 
QCD 

In this extra section, contained only in the 
arXiv version of this paper, we provide some pre- 
dictions from the AdS/QCD approach. 

In the holographic model of Brodsky and de 
Tcramond (see, e.g., [iQ]), the pion DA has the 
form 

^ V^(T^ (5) 
with the normalization 

= 5(3/2,3/2) = ^ (6) 

and moments given by 

hoi ^ 1 i?(3/2,(2n+l)/2) 

4 5(3/2,3/2) ' ^ ' 

where B(x,y) is the Euler Beta function. Taking 
the first derivative of (p{x)^°^, ~ one real- 
izes that this DA has its endpoints strongly en- 
hanced, even relative to the asymptotic DA, let 
alone to the BMS DA (cf. Fig. [fj). Note that 
also the twist-4 DA contribution, extracted from 



the holographic model by Agaev and Gomshi No- 
bary |41| . shows endpoint enhancement, although 
this effect is milder compared to the twist-2 DA 
case. It is interesting to recall at this point that 
Eq. ([5]) was considered long ago by Mikhailov and 
Radyushkin in an attempt to reconstruct the pion 
DA from its first few moments within the context 
of QCD sum rules with nonlocal condensates [6]. 

The values of the moments {£,'^)tt and ('f'*)^ of 
the holographic model have been inserted in Ta- 
bles [Tj and [2l respectively, (last entry in each of 
these tables) after NLO evolution from their ini- 
tial values (C^)^ = 1/4 and {^'^)^ = 1/8 to the lat- 
tice scale ^i^^ (with associated Gegenbauer coeffi- 
cients 4°'(Mfat) = 0.107 and a5°VL) = 0.038). 

The holographic model is also included in Fig. 
[2] in comparison with other theoretical models 
and confronted with the CLEO data [11], in par- 
allel with the two most recent and precise lat- 
tice calculations |18|19j . As for the moments 
above, we have assumed that the holographic 
DA, given by Eq. ([5]) is normalized at the scale 
/i^ = 1 GeV^ with the Gegenbauer coefficients 
al^°\^i^ = 1 GeV^) = 7/48 and a\°\n^ = 
1 GeV^) = 11/192. We then performed a two- 
loop evolution to the scale A^cleo ^^'^ found 

«2°'(a*cleo)= 0.101, 4-'(Mcleo) = 0.03 . (8) 

The position of the holographic model is marked 
in Fig. [5] by a red bullet and turns out to be (cf. 
Fig- E]) just inside the border of the 2a ellipse 
of the CLEO data, being also inside the predic- 
tions' range of the QCDSF/UKQGD Collabora- 
tion, but outside the limits on 02 determined by 
the UKQCD-RBC07 Collaboration. 

Predictions for the spacelike (and timelike) 
pion form factor, derived from AdS/QCD, can 
be found, for example, in [40|41|42|43j : they are 
not addressed here. 
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